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Abstract 

A four-parameter potential is analyzed, which contains the three dimensional harmonic 
oscillator as a special case. This potential is exactly solvable and retains several characteris- 
tics of the harmonic oscillator, and also of the Coulomb problem. The possibility of similar 
generalizations of other potentials is also pointed out. 


1 Introduction 

Searching for exact solutions of the Schrodinger equation has been an interesting challenge since 
the early period of quantum mechanics. This classic area has gained new momentum from the 
recent introduction of supersymmetric quantum mechanics (SUSYQM) [1], which relates pairs of 
essentially isospectral potentials to each other by means of (super)algebraic manipulations. (See, 
for example [2] for a recent review on SUSYQM, and [3] and references for its relation to other 
methods of analyzing isospectral potentials.) This new approach helped to view old problems 
from a new angle, and allowed unified, systematic treatment of previously unrelated results. New 
solutions of the Schrodinger equation have been described and classified, together with already 
known ones. The most well- known potentials have been shown to have the property of shape 
invariance [4], a concept introduced in SUSYQM. Much less is known, however, about the more 
general Natanzon potentials [5], which are, in principle, solvable, nevertheless their practical use 
is hindered by their complicated mathematical structure. The techniques inspired by SUSYQM 
allow a straightforward generalization of the simplest shape-invariant potentials, while avoiding 
most of the mathematical complications characterizing the general Natanzon potentials. 

Here I discuss a potential which can be considered the simultaneous generalization of the three 
dimensional harmonic oscillator and Coulomb potentials: these two shape-invariant potentials 
can be obtained from it by tuning one of its four parameters. Its ('oulomb limit lias already 
been described [6], and here I discuss its connection with the harmonic oscillator. I 11 contrast 
with other anharmonic oscillators, this potential converges to a finite value in the v * 00 limit. 
It also inherited several characteristics from its two “parent potentials’', which may enable its 
applications to physical problems, where deviations from these two fundamental potentials are 
relevant. 

I 11 Section 2. 1 give a brief account of a simple procedure which can be used to derive exactly 
solvable potentials. Section 3. contains the main results of this contribution, while in Section 4. 
a summary is given and directions towards hither investigations arc pointed out . 
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2 Transformations of the Schrodinger Equation 

Here I describe an old method of solving the Schrodinger equation to demonstrate how a wide 
range of solvable potentials can be derived in a relatively straightforward way. Originally this 
procedure was used [7] to derive only some well-known potentials, but it can be shown that the 
general Natanzon potentials can also be derived from it. This procedure was also connected to 
the formalism of SUSYQM [8]. 

The solutions of the one-dimensional Schrodinger equation (with h = 2m = 1) 


d 2 vj/ 

dr 2 


+ (E 
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( 1 ) 


are generally written as 

*(*) = f{*)F[g(x)), 

where F{g) is a special function which satisfies a second-order differential equation 


( 2 ) 


^ + Q(9)^ + K{<l) F (<l) = Q- (3) 

Here Q(g) and R(g) are well-known for any specified special function F(g), while /(. t) and g(x ) 
are some functions to be determined. Substituting (2) in (1) and comparing the results with (3) 
we arrive at the following expression [8] after some straightforward algebra: 
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Eq. (4) relates the only undetermined function //(./■) to the dilference of t he energy E and the 
potential V(x). Observing that the energy term E on the left -handside of Eq. (4) represents 
a constant, the authors of Ref. [7] equated certain terms of the right- -handside with a constant 
to account for it. This results simple differential equations for g(x). The authors in Ref. [7] 
applied this method to the hypergeometric and confluent hypergeometric function and obtained 
the solutions of some simple potentials. 

Considering the particular example of the confluent hypergeometric function F(— n, f)\ g(x)) 
and introducing the simple g{x) = ph(x) substitution we get 


E n - V(x) 


h m {x) 3 
2 h'(x) ~ 4 
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Identifying one of the last three terms on the right -hand side of Eq. (5) with a constant, the three 
shape-invariant potentials of the confluent hypergeometric case, the three-dimensional harmonic 
oscillator, the Coulomb problem and the Morse potential, are recovered. These potentials appear 
in the radial Schrodinger equation, therefore in what follows I shall replace x with r. 
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3 Generalization of the Harmonic Oscillator Potential 


A straightforward way of generalizing the simplest possible solvable potentials to more general 
ones is identifying combinations of several terms on the right-hand side of (5) with a constant. 
This procedure recovers the Natanzon confluent potentials [5], the solutions of which contain 
confluent hypergeometric functions. The most general six-parameter version of these potentials 
can be obtained by considering the combination of all three terms on the right-hand side of (5), 
which explicitly contain parameters, however, the technical difficulties increase considerably in 
this case. The problem remains relatively easy to handle if we take the combination of two such 
terms only. Considering the differential equation 

< v( *»’( i+ *£j)= c (6) 


corresponds to “mixing” the harmonic oscillator and the Coulomb potentials: 0 — > 0 recovers the 
latter one [6], while 6 — ► oo combined with C — CO yields the former one. (See Eq. (5)). The 
Coulomb limit has been discussed in detail in Ref. [6], and here we focus on the harmonic oscillator 
limit. The potential described here and in Ref. [6] is essentially the same for any finite' value of 0 , 
nevertheless, it is more convenient to use different notations when we discuss its connection to the 
two limiting case. In order to make the formalism of the two limits compatible with each other, 
here we follow the notations of Ref. [6] as closely as possible. 

As described in [6], the differential equation (6) can be solved explicitly for the inverse r(h) 
function only: 


r 




l # 1 / 2 tanh " 1 



+ 



(?) 


This function, of course, can be used to determine h(r) as well to any desired accuracy. 



FIG. 1. The h(r) function defined by Eq. (7), displayed for 0 = 0 . 1 , 1 , 10 , 100, oo 
and C = 1. (Curves lying higher correspond to higher value of 0.) 
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We have plotted h(v) in Fig. 1. for several values of the parameter 9. As discussed in Ref. [6], 
h(r) can be approximated by Cr 2 / 4 near the origin, and asymptotically follows h(r) — ► (C0) 1/2 r 
in the r — > oo limit, which correspond to the h(r) functions characterizing the harmonic oscillator 
and Coulomb problems, respectively. (See e.g. Ref. [8].) The range of the transition between 
these two regions is governed by the 9 parameter: it moves towards larger values of r as 9 increases 
(see Fig. 1.), and disappears completely in the 6 — + oo (harmonic, oscillator) limit. 

Substituting h(r) into (5) and removing the n-dependence from the potential terms by intro- 
ducing the constant 

__ j „ / n\ 

( 8 ) 
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(which amounts to a specific choice of p = />„,) we arrive at the following potential 
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These formulas differ from the corresponding ones in Ref. [6] only in a shift of the energy scale 
and in the usage of slightly different parameters. The changes reflect the difference between the 
Coulomb and harmonic oscillator limits of the general problem containing both potentials as a 
special case. These differences, however, do not essentially influence the form of the wavefunctions: 


CM**? l r (» + /i) \ 1/2 

r(^) \ n -((/^ + 2n)0 _] + pn)J 

x ( 1 + /i(r)/0) 1/4 (/r(r))^ zi exp F{-n, /?; p n k(r)). (11) 


(Here and in Eq. (10) n denotes the number of nodes in the radial wavefunction.) 

As we can expect from (6), these formulas reduce to the corresponding ones for the harmonic 
oscillator in the 9 — > oo limit, if we introduce the notation u> = CD'/ 2 and 1 = 0 — 3/2. In 
particular, the two last terms in (0) vanish and the first and second terms transform into the 
harmonic oscillator and centrifugal terms, respectively. We have displayed V (?•) and the position of 
some of the lowest- lying energy eigenvalues in Fig. 2. for some values of parameter 9. As it can be 
seen there, the oscillator character of the potential strenghtens with increasing 9. V{i') is oscillator- 
like near the origin, and approximates the Coulomb potential (with Ze 2 = C 1 / 2 DO 3 / 2 ) for large r. 
The domain of oscillator-like behavior expands with increasing 9: this is related to the structure 
of h(r ) discussed previously. (See also Fig. 1.) Also, the energy spectrum is oscillator-likejor 
small values of «, and Coulomb-like for large n: E n converges to E n ^^ = V(r — > oo) = CDO. 
See Ref. [6] for a more detailed description of V'(r) in terms of powers of r. 
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FIG. 2. Potential V(r) in Eq. (9) displayed together with the lowest-lying energy 
levels for 6 = 1, 10, 100 and _oo. The other parameters are C = 1, D = 5 and 0 = 1.5 
in all cases. V{v — * ex)) = C DO in each case. 

Similarly to the Coulomb limit discussed in [6], this potential can be rewritten into the sum 
of a central, centrifugal and /-dependent part: 

/(/ + 1 ) 

V ( r ) = V'o(r) + V/(r) + » 


( 12 ) 
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The definition of / is, however, different in the two limits: / = ft / 2 — 1 for the Coulomb [6], and 
l — ft — 3/2 for the oscillator limit. (Fig. 2. displays potentials with / = 0 only.) Also in contrast 
with the Coulomb limit, v(r) in Ecp (14) does not vanish for large values of ?■, rather it goes to 
the value —3/4. This, again, is the consequence of the asymptotical Coulomb-like character of 
V(r). 

It is remarkable, that E n depends on the combination 2 n + ft only (i.e. on 2n + / + 3/2 
in the oscillator limit), therefore the generalized harmonic oscillator potential has a degeneracy 
pattern similar to that of the harmonic oscillator. In other words, the terms representing the 
anharmonicity do not remove the degeneracy of the energy levels. 

This generalization of the harmonic oscillator potential could be applied to physical problems, 
where an attractive Coulomb potential is distorted by an oscillator-like potential component for 
small values of r. I his is the cast', for example, for a finite, homogenous, spherical charge distri- 
bution, but in that case the resulting potential can strictly be separated into two domains, where 
it exactly follows r 2 like and v 1 like behavour. The potential discussed here can be considered a 
deviation from this simple model problem. An example for a similar situation is discussed in Ref. 
[9] in connection with a potential experienced by electrons in certain crystal environments. 

Finally, there are some other potentials occupying a similar intermediate position between the 
simple shape invariant potentials and the general Natanzon potentials. Some of these, like the 
Woods-Saxon [10] and Ginocchio [11] potentials have been found earlier, while some others, the 
‘PHI [12] potential and those in Reis. [13,14,15] have been identified only recently, mainly in 
SIISYQM related studies. See Ref. [6] for more details. 


4 Summary 

Here I have analyzed a four- parameter potential, which contains both the harmonic oscillator and 
the Coulomb potential as special cases. I have interpreted this potential as the generalization of 
the harmonic oscillator potential, and have established that it is a special admixture of a long- 
range attractive Coulomb term, and an oscillator— like term near the origin. This is also reflected 
in the structure of the energy spectrum. 

Exact analytical solution of the radial Schrodinger equation can be obtained for any partial 
wave, however, an angular-momentum-dependent term appears for / / 0. A remarkable finding 
is that the anharmonicity appearing in the general form of the potential does not remove the 
degeneracy of the energy levels. 

Similar generalizations of the harmonic oscillator and other well-known potentials are also 
possible by considering further simple differential equations similar to that in Eq. (6). These 
subclasses of the Natanzon potentials seem to be suitable for applications, because they have 


118 


more flexible shape than the simplest solvable potentials, but may still remain relatively simple 
to handle mathematically. 
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